Charge inversion at minute electrolyte concentrations 
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Anionic DMPA monolayers spread on LaCl3 solutions reveal strong cation adsorption and a sharp 
transition to surface overcharging at unexpectedly low bulk salt concentrations. We determine the 
surface accumulation of La 3+ with anomalous x-ray reflectivity and find that La 3+ compensates the 
lipid surface charge by forming a Stern layer with ~ 1 La 3+ ion per 3 lipids below a critical bulk 
concentration, c t » 500 nM. Above Ct, the surface concentration of La 3+ increases to a saturation 
level with « 1 La 3+ per lipid, thus implying that the total electric charge of the La 3+ exceeds 
the surface charge. This overcharge is observed at ~ 4 orders of magnitude lower concentration 
than predicted in ion-ion correlation theories. We suggest that transverse electrostatic correlations 
between mobile ions and surface charges (interfacial Bjerrum pairing) may account for the charge 
inversion observed in this dilute regime. 

PACS numbers: 73.30.+y, 82.45.Mp 



Counterion screening of charged interfaces in elec- 
trolytic solutions is key for understanding a broad range 
of phenomena in molecular biology, colloid and polymer 
science, or microfluidics. Directly or indirectly, charges at 
cell membrane surfaces control functions and conforma- 
tions of nearby molecules that may be involved in inter- 
and intracellular transport processes, cell-cell recognition 
and biomimetic mineralization processes. Despite exten- 
sive experimental and theoretical work performed over 
more than a century, problems regarding the nature of 
ion correlations, the role of the hydration sheath in phys- 
ical processes and the structure of water near charged 
interfaces remain the subject of intense theoretical and 



experimental work [1|, |2j, 



3j. 



The phenomenology of the electrostatics in soft media 
is rich and complex. Effects such as charge inversion, 
where charges at an interface attract counterions in ex- 
cess of their own nominal charge density have been re- 
ported in the literature 0, 0, II IE El. Current theories 
assume that charge inversion results from the free en- 
ergy gain brought about by ion-ion correlations [l| and 
predict that measurable charge inversion occurs when 
q/ (iter A 2 ^) ~ 1, where q = Ze is the counter-ion charge 
and a is the surface charge density. Xd = (87r/s/) _1 / 2 
is the Debye length of an electrolytic solution (Bjerrum 
length, Ib = e 2 /(efceT)) with ionic strength / |7J. Ap- 
plying this reasoning to Langmuir films of a charged lipid 
such as phosphatidic acid, one expects for a typical area 
per amphiphile, An p id(— e/er) « 40 A 2 , spread on LaCi3 
solution {Z = 3), that the salt concentration has to be at 
least in the 10 mM range for charge inversion to become 
significant. Besides correlations of the mobile ions, there 
are other effects that may contribute to the phenomenon. 
For example, it has recently been suggested that trans- 
verse correlations, i.e., correlations between interfacial 
charges and mobile ions, may also play a role in the gen- 
eration of charge inversion [8|, |9j . 



In this paper, we demonstrate charge inversion at 
ultra-low (/xM or less) LaCl3 electrolyte concentrations 
with surface-sensitive, resonant X-ray scattering at the 
La absorption edge and discuss its relationship to 
proton transfer and release. These results extend previ- 
ous studies of BaCL [ijJIlll and Csl solutions [H[l3|. 

Surface monolayers of DMPA (1,2-dimyristoyl-sn- 
glycero-3-phosphatidic acid) (Avanti Polar Lipids), were 
spread from chloroform/methanol (3:1; both from Merck, 
p. a. grade) on aqueous solutions with different LaCl3 
(Sigma, purity: > 99%) solutions prepared from ul- 
trapure water (NANOpure, Barnstead) in an enclosed 
Teflon (PTFE) Langmuir trough. For the x-ray scatter- 
ing work, the air space above the trough was continuously 
purged with water-saturated helium. Both isotherm and 
scattering studies were performed at a temperature of 
21°C. To minimize ion contamination, the subphases 
were handled and transferred to the Langmuir film bal- 
ance in bottles made from Teflon (Fisher Scientific). Af- 
ter solvent evaporation, the monolayers were compressed 
at a rate of ~ 1 A 2 /(molecule min.) to record isotherms 
or measure reflectivities at various surface pressures be- 
tween 7r = 10 and 50mN/m. At bulk ion concentrations 
in the /zM range, monolayers become progressively un- 
stable to compression. At a collapse pressure, 7r c , the 
compressibility of the surface film increases, consistent 
with a gradual transformation of the monolayer locally 
into bilayer or trilayer structures. n c decreases with in- 
creasing salt concentration, which makes measurements 
at high 7T difficult for c^ h > 1 fM. We restrict the 
following discussion to the low surface pressure regime, 
tt = 15...30mN/m. 

X-ray reflectivity (XR) measurements were conducted 
on the Ames Laboratory horizontal surface diffractome- 
ter at beamline 6-ID-B of the Advanced Photon Source 
(APS). The highly monochromatic beam (energy resolu- 
tion, AE ~ 1 eV) , selected by an initial Si double crys- 
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tal monochromator, is deflected onto the liquid surface 
at a specified incidence angle by a secondary Ge(lll) 
monochromator located at the difffactometer The 
x-ray energy E was calibrated by measuring the absorp- 
tion spectrum of a dilute LaCl3 solution. Off-resonance 
spectra were measured at 8.0 keV, and "on-resonance" 
refers to the La L]jj edge at E res = 5.486 keV. To re- 
duce radiation damage, the Langmuir trough was period- 
ically translated across the beam. To monitor sensitively 
for damage, in particular after prolonged exposure of the 
sample at high momentum transfer, Q z , measurements 
were routinely repeated - without translating the sam- 
ple - across the Q z regions of the sharp cusps in the 
interference minima. 

To extract scattering length density (SLD) distribu- 
tions across the interface from the XR data, parame- 
terized profiles p(z) = p'(z) + ip"{z) were constructed 
in which the real and imaginary parts of p describe 
the electron density (ED) and the absorption density 
(AD), respectively, along the surface normal^z. We used 
both a modified 'box' model approach [13|, [l5| and the 
quasi-molecular Volume-Restricted Distribution Function 
(VRDF) approach [U 13 for this reconstruction to ob- 
tain stable results and to check for consistency. In the 
box model, p(z) is described by a sum of error functions 
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where N is the total number of slabs j that describe 
the monolayer. Zj and o~j denote the position and r.m.s. 
roughness of the jth interface, respectively, between the 
slabs, and j = N + 1 indicates the bulk subphase. Alter- 
natively, p(z) is constructed from a sequence of thermally 
broadened distributions, modeled as Gaussian functions, 
that represent parsed subfragments of the lipid structure 
and bound cations [Til In both approaches, volume 
filling is implemented by accounting for the subfragment 
volumes as derived from molecular dynamics simulations 
[I]) . For either model, p"{z) is only significant in spec- 
tra measured at E res and correlates with the physical 
distribution of the cations across the interface, i.e., it is 
negligible where the cation density is at bulk density or 
below. Reflectivities are calculated from p{z) by apply- 
ing a recursive method [l~8j to the discretized profiles. 

Both / Ap'(z) dz = J (p' off _ res {z) - p' res {z)j dz, and - 

more directly - J p"(z) dz provide a direct and model- 
independent measure of the number of La 3+ adsorbed at 
the interface: 



/ Ap' dz 

Az e ff ■ Alipid 



(2) 



where Az e g ~ 19.0 is the difference between the effective 
numbers of electrons for La 3+ at E jj- res and E res . 

Figure ^ shows isotherms, tt vs. Au P id, of DMPA 
monolayers at various LaCl3 concentrations of the bulk 
subphase. On the salt-free subphase, a sloped plateau in- 
dicates coexistence of a liquid and a hexatic lipid phase. 
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FIG. 1: DMPA monolayer isotherms for different La 3+ sub- 
phase concentrations 
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FIG. 2: Fresnel-normalized x-ray reflectivities, R/Rf, of 
DMPA monolayers at E = 8.0 keV (off-resonance) and 
5.486 keV (La Liu resonance) on aqueous LaCl3 solutions. 
For clarity, subsequent pairs of data sets are offset by a factor 
of 100 each. Solid lines derive from model fits. 



The height of this plateau depends sensitively on salt 
concentration. As shown earlier, micromolar concen- 
trations of Ca 2+ condense the monolayer to the extent 
that this plateau disappears 0] . Figure Q shows that 
trivalent cations are even more effective in monolayer 
condensation and the plateau has already been lost at 

„LaCl 3 



hulk 



> 10 nanoM, consistent with strong electrostatic 
interactions between La 3+ and the DMPA - headgroups. 
It has been shown with fluorescence microscopy that 
Ca 2+ condenses DMPA - from the two-dimensional (2D) 
gas phase into hexatic domains 0|- We expect that 
binding of La 3+ to DMPA - has a similar effect already 
at exceedingly low salt concentrations. 

Figure [21 shows exemplary, Fresnel-normalized off- 
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FIG. 3: VRDF electron density profiles derived from fits to 
the experimental data. The off-resonance and on-resonancc 
data shown in Fig.|5|were co-refined with consistent sets of pa- 
rameters that deviated only in the effective electron numbers 
and absorption cross sections of the La 3+ component. Subse- 
quent ED pairs are offset by 0.2 e~/A 3 for clarity. Similar ED 
profiles were derived from a modified box models as described 
in the text. The number density of La 3+ at the interface was 
subsequently obtained from such models via Eq. 



resonance/on-resonance reflectivity pairs of DMPA 
monolayers at 7r = 15mN/m on subphases that contain 
100, 500, and 1 fiM LaCi3. A high-energy reflectivity of 
DMPA on pure water is also shown. Qualitatively, the 
increase in intensity with increasing LaCl3 bulk concen- 
tration indicates ED increases due to the accumulation 
of La 3+ at the acidic interfaces. Quantitatively, the dif- 
ference between off-resonance and on-resonance EDs is 
proportional to the La 3+ surface concentration (Eq. [3J . 
The data fall in two categories: Between c£°^ 3 = 
and 500 nM LaCi3, both off-resonance and resonance re- 
flectivities are rather similar. Above 500 nM LaCkj, on 
the other hand, the off-resonance reflectivities are signif- 
icantly higher at all Q z than those at resonance. With- 
out any model interpretation, this qualitative difference 
shows conclusively that the condensation of La 3+ at the 
monolayer is significantly larger at 1 fjM LaCl3 than at 
500 nM and below. This is corroborated and quantified 
by a detailed data analysis (see Fig.^J, using the modified 
box model and the VRDF. We also measured more com- 
plete, systematic sets of off-resonance reflectivities over 
a wider range of subphase salt concentrations, between 
100 nM and 10 fjM, and at higher surface pressures. As 
shown in Fig.^J these results are consistent with the more 
stringent anomalous reflectivity data. 

Figure |3| shows SLD profiles p(z) obtained from a 
co-refinement of data pairs shown in Fig. [21 While 
p'(z) in the lipid headgroup region is larger for the off- 
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FIG. 4: Number of adsorbed La 3+ per DMPA in surface 
monolayers as a function of bulk LaCi3 concentration. Red 
plot symbols indicate results at 15 mN/m derived from the 
anomalous reflectivity data pairs shown in Fig. |5]and ED pro- 
files similar to those shown in Fig. [3] The results indicated by 
blue symbols were derived from independent data measured 
at 8.0 keV alone. 



resonance data than for the on-resonance data on all 
LaCl3-containing subphases, the difference is much larger 
for 1 fjM LaCl3. The VRDF modeling suggests that the 
differences in p', to be interpreted as the distribution of 
cations at the interface, are confined to the lipid head- 
group region. Equation [3 serves to determine n La 3+, the 
number of La 3+ per DMPA headgroup in the monolayer. 
A compilation of this data collected at various surface 
pressures as a function of ' 3 , is given in Fig. 0J The 
plot shows clearly three distinct regimes which consist of 
two plateaus at low and high c^°^' 3 , separated by a sharp 
transition. In regime I, between 100 nM and 200 nM 
LaCl3, n ha 3+ is constant at 0.34 ± 0.1. It is followed by 
a remarkably sharp transition, regime II, that connects 
this regime with regime III, observed at high LaCl3 con- 
centration (> 1 /iM), where n La 3+ ~ 1. 

The results in regime I correspond to « 1 La 3+ bound 
to three DMPA at the interface. This stoichiometry fol- 
lows from the dissociation constants of the first and sec- 
ond protons to DMPA, pK a ,i = 2.1 and pK Q:2 = 7.1 [2l|. 
and the binding constant of La 3+ to a monovalent phos- 
phate group, estimated as Kb ~ 100 M _1 |9j], introduced 
into the PB theory yj, |20| • Regime II follows from the 
large affinity of La 3+ ions to divalent phosphate groups 
at the interface, K B ~ 10 5 - 10 6 M -1 9], which leads to 
deprotonation, PC^H - — > P0 4 , of the interfacial group 

at a concentration c^°^ + ~ 10~ 6 . The increase of n La 3+ 
from 1/3 to 2/3 is therefore not due to charge inversion, 
but rather to a doubling of interfacial charge density by 
surface deprotonation. Regime III with ?i.La 3 + ~ 1, on 
the other hand, clearly shows charge inversion that is 
not covered by established ion- ion correlation theories . 
The experimental observation of charge inversion in this 
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concentration range is thus compelling evidence for yet 
other effects of the interface A hypothetical Stern 
layer, that would consist of La 3+ bound to interfacial 
PO4 - in a 1:1 ratio (n La 3+ 1), leads inevitably to a 
diffuse distribution of co-ions (e.g., CI - ) which results in 
a free energy penalty (per amphiphile) of Q 

F « 2(ln(2V2A D /A GC ) - 1) « l8k B T, (3) 

where Xqc = e/(27r<TZs) is the Gouy-Chapman length. 
We envisage two scenarios, under which such a free en- 
ergy penalty could be compensated by other effects: (1) 
A "molten salt" state consisting La 3+ ions intercalated 
between the negative charges on the phosphate oxygens 
forming a highly correlated state 0- (2) Enhanced ad- 
sorption in the Stern layer via non-electrostatic interac- 
tions, for example, hydrogen bonding between the PO4 - 
oxygens (acceptor) and a La(OH) 2+ complex (donor). In 
either scenario, charge inversion is associated with some 
instability of the Langmuir monolayer, as evidenced by 
slow collapse at intermediate surface pressures at LaCl3 
concentrations above 5 • 10 -5 M. 

In this paper we presented experimental evidence for 
charge inversion at minute electrolyte concentrations 
(0.5 fiM LaChj) in contact with acidic DMPA monolay- 
ers. When compared with experimental work investigat- 
ing forces between charged plates |(|, charge inversion 
in our experiment is observed at « 2 orders of magni- 
tude lower and in the absence of any monovalent salt. 
More fundamentally, our results provide strong indica- 
tion that correlations among mobile ions and interfacial 
charges play a prominent role in the generation of charge 



inversion. We interpret the concentration dependence of 
counter-ion accumulation beyond regime (I) in Fig. 0] 
as follows: Charge inversion is initiated by the deproto- 
nation of the secondary oxygen in the phosphate group 
(regime II) . The high level of charge inversion at fjM bulk 
LaCl3 (regime III) requires additional effects of ion reor- 
ganization, and we suggested two tentative scenarios. 

The results presented in this study corroborate the no- 
tion that charge inversion is a generic effect for multiva- 
lent ions next to highly charged surfaces and quantify the 
role of multivalent ions in favoring surface deprotonation. 
Furthermore, however, these results are likely of univer- 
sal relevance for the physics of highly charged colloidal or 
biological systems in the presence of multivalent ions, as 
evidenced from previous studies HG3- Finally, it is also 
clear that still more experimental and theoretical work 
is needed to fully resolve the fundamental issues involved 
in the electrostatics of highly charged interfaces. 
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